Non-Abelian vortices in the emergent U(2) gauge theory of the Hubbard model 
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By the spin-fermion formula, the Hubbard model on the honeycomb lattice is represented by a 
U(2) gauge theory in the mean field method, non-Abelian vortex solutions are constructed based 
on this theory. The quantization condition shows that the magnetic flux quanta are half-integer. 
There are 2k bosonic zero modes for k winding vortices. For the fermions, there are 2 zero energy 
states (ZESs) corresponding to the single elementary vortex. In the vortex core and on the edge, 
the system are in the semi-metal phase with a spin gap and in the insulator phase with Neel order 
phase, and can be mapped to the superconductor in class A and CI, respectively. 

PACS numbers: 71.10.Fd, 64.60.-i 



In 2+1 dimensional systems, electrons moving on the 
lattice can be well described by the Hubbard model. It 
is a long-standing research topic due to the richness of 
its phase diagram and its possible explanation of high-T c 
superconducting. The Hubbard model on the honeycomb 
lattice is also suitable to study novel materials like the 
graphene, which has a great potential for applications. 
Previously, Hou, Chamon and Mudry (HCM) discussed 
the Abelian vortex in graphene-like structures^, where 
the vortex configuration of the order parameter is caused 
by the distortion on the Kekule texture. Herbut used the 
Hubbard model to consider vortices in the background 
Neel order, and found two orthogonal ZESs of fermions at 
the vortex core^. Nevertheless, the vortex configurations 
in the above cases are of Abelian type. 

In recent years, vortices of non-Abelian type have 
been discovered in M = 2 supersymmetric QCD 
(SQCD)2r— , which is indeed a breakthrough since the 
Abelian Abrikosov-Nielsen-Olcsen (ANO) vortex has 
been constructed-,. The similarity of the emergent gauge 
theory (EGT) of the Hubbard model and the bosonic 
truncation of SQCD motivates us to construct the non- 
Abelian vortices in this theory. In this work, we present 
a complete investigation of the EGT of the Hubbard 
model in the low energy limiti. By the mean-field the- 
ory method, we construct the non-Abelian vortex solu- 
tions in the theory. Non-Abelian vortices have non-trivial 
orientational zero-modes, which induce net interactions 
between multi-vortices^. We hope that the non-Abelian 
vortices enable us to understand novel phenomenons in 
condensed matter physics, especially the unconventional 
superconductor. 

The honeycomb lattice has a special property that the 
occupied and empty states meet at two Fermi points in 
momentum space. Due to the bipartite property of the 
Honeycomb lattice, one can introduce three Pauli ma- 
trices a a , r a , and p a to act on the spin, the sublattice 
(pseudospin) and the valley spaces, respectively. The 
standard Hamiltonian of the Hubbard model is written 



as H t + Hu, where H t describes the hopping of electrons, 
and Hjj describes the repulsive interaction of electrons. 
Expanding H t near the Fermi points, one obtains 9 

Ht = i^J ^kClj^ir^ + r^ky)^^ (1) 

where t is the hopping constant. The index s = (t,i) 
denotes the spin, p,q = (A,B) denote the sublattice, 
and i,j — (1,2) denote the valley. With Hjj term, the 
Hubbard model has a very rich phase diagram. 

The emergent gauge theory kicks in an elegant way 
to unify phases of the Hubbard model. Hermele con- 
structed an SU(2) gauge theory for the Hubbard model 
on the honeycomb lattice, where the electron operators 
were written in the slave-rotor formulation, and an SU(2) 
algebra spin liquid (ASL) is found on the insulator side of 
Mott transition 10 . Sachdev et al. constructed an SU(2) 
gauge theory of the Hubbard model in a different way^^, 
which will be presented in the following. By making use 
of the spin-fermion formula, the electron operators can 
be decomposed as 

where R denotes the spin density wave (SDW) order. ip± 
are the spinless fermion operators, the indices ± measure 
the spin-projection along the local spin reference axisS. 
The R fields can also be interpreted as the " auxiliary bo- 
son" in the formulation of Kotliar and Ruckenstein (KR) , 
which represents four boson operators corresponding to 
four atomic states per site 11 . In Ref 0, a Higgs-like po- 
tential has been constructed for R, which leads to sym- 
metry breaking and is necessary for vortex solutions. In 
this manner, R is generalized to an invertible complex 
2x2 matrix, which represents a generic spin operator for 
arbitrary filling^. 

Let us show how the hidden gauge symmetry emerges 
from the spin-fermion formula. First, the R matrix can 
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TABLE I. The transformation properties of the fields. 
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(R) = 


Semi-metal 


Insulator with VBS order 


(R)^0 


Scmi-mctal with a spin gap 


Insulator with Noel order 



TABLE II. The phase diagram of the system with U(2) gauge 
theory^, in which VBS is the abbreviation of valence bond 
solid. 



be rotated by an SU(2) global spin rotation from left 



c ^VC 



(3) 



As this action is global, the spin rotation will be re- 
garded as the "flavor" symmetry. Secondly, there is a 
local gauge redundancy, the theory is invariant under the 
gauge transformations 



C^C 



(4) 



where U is an element of a gauge group. Sachdev et al. 
consider U € SU{2) , here we argue that the system is 
invariant under C7(2), and the theory is similar to the 
Weinberg-Salam (WS) electro-weak theory, but limited 
to 2 + 1 dimensions. This extra U{1) gauge group is the 
key ingredient for the non-Abelian vortices, because the 
homotopy group is non-trivial for t/(l) x SU(2), i.e., 



(5) 



Instead, when the SU{2) gauge group is considered, only 
TL-i vortex can be constructed, which is genuinely of 
Abelian typei^. 

The Lagrangian of the emergent U (2) gauge theory of 
the Hubbard model reads^ 

C =*7 M (5 AI + io*A%)V - A$ a *pV a * 

+ ~ [(d^ a - 2e abc A b ^ c f] + S ($ a ) 2 + u ($ a ) 4 

+ Tr[2? A1 i?(2? Al i?) t ] +sTr(R jf R) + u[Tr(R^R)] 2 , 

(6) 

where the covariant derivative is defined to be T>^R = 
dfj,R — iA^Ra a . The index a runs from to 3, mean- 
while a, b take value as a, b = 1, • • • , 3 in the following. 
"0" stands for the Abelian {7(1) part of the gauge group, 
and a, b stand for the non-Abelian SU(2) part. The rep- 
resentation of the algebra is that <t° = I2, <J a is the Pauli 
matrix, which follows that Tr[a a a p ] = 2S al3 . Hence $ a 
transforms as the triplet of gauge SU(2), p z is the third 
component of the Pauli matrices which acts on the val- 
ley space, see Eq.([T]). The theory has a set of symmetries 
U(l) g ® SU(2) g <g> SU(2) S , the transformation properties 
of the matter contents $ a and R are listed in Table fl] 
The s, u and tildes of them are the coupling con- 
stants which tune the couplings of R, $ a . Four differ- 
ent phases are unified in the Lagrangian of Eq.([B]). They 
are classified by whether the vacuum expectation val- 
ues (VEV) of the R and $ are zero or noti, see Table 



III1 When the adjoint scalar <£> a is condensed at a energy 
scale v\ = —s/2u, the gauge symmetry breaking patten is 
SU(2) — > £7(1). the Lagrangian in Eq.© reduces to the 
f7(l) gauge theory, where monopole events (or hedgehog) 
occur 14 . When $" acquires a VEV, the Yukawa coupling 
term — \^ a ^> p z a a ^ will become a mass term for the 
effective mass is = \f2\v\. We decouple $ a by a 
hierarchical symmetry breaking, i.e. V\ 3> t>2, where 1)2 
is the VEV of R. The spectrum of <& a will not appear 
below the low energy scale i>2- 

R potential in Eq.© is generic, we choose a special 
potential of R in the following, which is written as 



V(R) = wTr 



1 



~-vii 



(J) 



where v\ = —s/u. Up to gauge rotations, the VEV of R is 
(R) = W2/v / 2l2, which is called the "color-flavor" locked 
(CFL) phase. In the Higgs phase, the masses of the two 
flavors (spin indices) are the same. In the CFL phase, the 
U(l) x SU{2) gauge symmetry is broken. Nevertheless, 
there is a diagonal invariant rotation, that is 



V f (R)Ut = (R), V f , L7 c t e SU(2). 



(8) 



This means the spin rotation of R can be equally done 
by the gauge rotation in the vacuum. So a combined 
color-flavor SU(2) c+ f global symmetry remains. 

With all the criteria above, at an energy scale below 
vi, the effect action can be written as 



(9) 



S*=J d 6 x{^V^ + -m&tf/jV*}, (10) 
Sn= J d 3 xTr{^\F^\ 2 + \V^R\ 2 + V(R)}, (11) 

where F^ = i[-D^, D v \, and g is the coupling constant 
of gauge fields. The Yang-Mills term is added by the 
consideration of relativistic dynamics. The static vortex 
configuration of the system is the minimal energy bound. 

The action in Eq. (fTTj) is the non-Abelian generalization 
of the Abelian Higgs model. Setting f3 = 2u/g 2 , which is 
the Landau-Ginzburg parameter classifying the types of 
superconductor. We will work in the BPS limit j3 = 1, 
which means |w| = g 2 /2. Performing the Bogomol'nyi 
completion, one obtains the BPS equations of the system 



F12 



VR =0, 



V(i2tR-fla) =0, 



(12) 
(13) 
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where 21? = T>i + iD 2 and z — x\ + ix 2 is the standard 
complex coordinate. The energy bound of the system is 



£ > - J d 2 xTT(F 12 vl) = -2nvlk, k e %. 



(14) 



The generators of SU(2) are traceless, so it does not 
contribute to the bound energy. Meanwhile, the £7(1) 
part does contribute, and the minimal energy bound is 
the quantized magnetic flux. Notice that our winding 
number is k/2, which is a half-integer— The normaliza- 
tion condition, Tr(cr°) = 2, cancels the "2" factor on the 
denominator of the winding number. This can explain 
naturally the half-integer magnetic flux quanta found in 
certain superconductor materials. 

Choosing the minimal winding, the ANO-like Ansatz 
can be embedded in the upper left-hand corner in R , i.e, 



R = 



e* e fi(r) 

MO, 







(15) 



which is named the (1,0) vortex. Similarly, the wind- 
ing term can also be embedded in the lower right-hand 
corner, which is named the (0, 1) vortex. Both cases 
have degenerate energy. Here /i(r) and /2(f) are the 
profile functions, the boundary conditions for them are 
h (00) = / 2 (oo) = V2/V2, /i(0) = and r / 2 (O) = 0. 
Requiring T>iR — > when r —¥ 00, the Ansatz for the 
gauge field Ai is written as 

A i = [(! ~ 9x(r)) 1 2 + (1 - 52M) a 3 ] . (16) 

The profiles 51,2 satisfy the boundary conditions 
51,2(00) = and 51.2(0) = 1. One can see that it is the 
combination of U(l) and SU(2) which cancels the diver- 
gence of diR. Substituting the Ansatz into the Eq.(fT2|) 
and (TT3"j) . one has the BPS equations for the profiles 



d h , 

-T-fi 5i 

ar r 



52) -0, iA 5l - 5 2 (/ 1 2 + /. 
r ar 



^-h - f -{9i~ 52) = 0, ~g 2 - g 2 (fl /I) = 0. 
ar r r ar 

(17) 

These equations can be solved numerically. 

In Eg. (ITS)) , the configuration of the (1,0) vortex breaks 
the global SU(2) c+ f color-flavor symmetry down to 
U(l) c+ f € SU(2) c+ f , the Nambu-Goldstone-like modes 
arise, which are named as the internal orientational zero 
modes, expressed as^ 



CP 1 



SU(2) c+f 



c+f 



(18) 



Other solutions can be generated by " color-flavor" trans- 
formations, i.e., R U c+f RUl +f , Ai -s- [7 c+/ A i C/ c t +/ , 
where U c+ f parameterizes the (DP 1 . All solutions have 
the same energy and the same boundary condition up to 



a regular gauge transformation. Thus, we have a mod- 
uli space. One can use the index theorem to calculate 
the dimension of the moduli space, which is found to be 
1 = N c Nf k/2 = 2fc2. These bosonic zero modes stand 
for the remanent excitation of order parameter, and are 
massive modes in the bulk regime (where R is in vac- 
uum). Therefore, they are localized in the vortex corei^. 
In the low energy, the dynamics of the zero modes can be 
well described by the (DP 1 sigma model^. Our analy- 
sis above is semi-classical. When quantum mechanism is 
considered, the (DP 1 sigma model will develop two vacua 
corresponding to the (1,0) and (0,1) elementary vortices, 
respectively^. 

Next, we will consider how many ZESs for the 8 com- 
ponents Dirac fermions. For Abelian vortices, k ZESs of 
fermions in the vortex-fcrmion system correspond to the 
fc-vortex background^ 7 -. The HCM model discussed the 
Abelian vortex in Kekule order parameters, see Eq.(3) in 
RefQ. When the spin indices are summed, two ZESs are 
expected for 8 components Dirac fermions. Turning on 
the Hjj term, Herbut also found two orthogonal ZESs at 
the center of the vortex for the Hubbard model on the 
honeycomb lattice 2 -, but the vortex configuration is in the 
Neel order N a , which exists in the spin space. The coinci- 
dence of the two instances is due to the fact that the Neel 
order N = (iVi ,^2,0) in Hjj will turn into an Kekule-like 
order A, i.e., A = Ni + iN 2 , when the the Hamiltonian 
is expressed explicitly by using the conventional 4 com- 
ponent Dirac spinor (itf,, u a , v a , Vb) as m 

RefQ. For our 

non-Abelian vortices, one of the two elementary vortex 
configurations (1,0) and (0,1) survives after considering 
the quantum effects, and can be set as the background 
vortex. Essentially, the (1,0) or (0,1) vortex is of Abelian 
type. Since R can represent the Neel order by the relation 
N a T a = Rt 3 R\ analogous zero energy Dirac equations 
in Herbut 's case can be given for Hamiltonian in Eq.([T]). 
So we expect two ZESs for the (1,0) or (0,1) vortex in the 
extended Hubbard model, and 2k ZESs for the k winding 
vortices. 

It is still an open question how the bosonic orienta- 
tional zero modes interact with the fermions in the vor- 
tex core, a self-consistent way to treat this question is to 
draw on the Bogoliubov-de Gennes (B-dG) equation 18 . 
In the B-dG equations, the fermion zero modes depend 
non-trivially on the winding term, it can be inferred that 
the degree freedom of the bosonic zero modes will not 
affect that of the fermion zero modes. It is shown that 
the fermion zero modes is well-localized in the vortex 
cor o 17 ' 18 . In ReffTil. singlet and triplet states of the 
unbroken symmetry are found for the non- Abelian vor- 
tices in high density QCD, where the SU(3) c+ l+r sym- 
metry is broken to SU(2) c+ l + r x U{\) c+ l+r- In our 
model, only U(l) c+ f symmetry group exists after sym- 
metry breaking, one singlet state can be constructed. We 
leave the B-dG equations of the extended Hubbard model 
for future work. 

It is interesting to notice that the vortex core and 
the vortex edge are in quite different phases. Far away 
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from the vortex core, the gauge bosons A" obtain the 
mass ~ 2gv2- The magnetic flux F\% cx 1/r 2 asymp- 
totes to zero. The non-existence of magnetic fields pre- 
serves the time reversal symmetry (TRS). With the set- 
up ($ a ) 7^ 0, $ a will turn into a Neel order parameter. 
The system is in the insulator phase with the antifcr- 
romagnetism (Neel) order 7 . If ($") = 0, the system is 
in the semi-metal phase. In the CFL vacuum, the spin 
rotation symmetry (SRS) is respected. The insulator 
with both TRS and SRS is associated with supercon- 
ductor in class CI according to Altland-Zirnbauer (AZ) 
classificatio n 19 ' 20 . 

In the vortex core, the physics is rich. ($ a ) ^ 
means that the SU(2) gauge symmetry is broken to 
£7(1) £ SU(2)1. The half quanta of magnetic flux pass 
through the core, so the TRS is broken. Considering the 
(1,0) vortex as the background, the boundary conditions 
of the profiles lead to B = diag(0, ^(O)) when r — > 0, 
where /2(C)) = 77 is a normalized constant. This indicates 
that the £7(2) gauge and SU{2) spin symmetry are not 
completely restored in the vortex core. Evidently, the 
SU(2) SRS is not only broken along the a 3 axis. This 
can be mapped to the superconductor in class A2i. A 
spin discrepancy occurs, i.e., Cf = 0, C± = Jpp-. The 
spin up fermions are suppressed, one can expect that 
the magnetism occurs in the core. If we take the (0,1) 
vortex, the magnetism will inverse the direction. In the 
presence of one pair of them, we will have a magnetic 
dipole in the system. Recall that our gauge symmetry is 
U{2) , an £7(1) x £7(1) gauge theory is not broken. With 
two spin flavors, the system becomes an £7(1) x £7(1) 
two-gap superconductor 22 . ^ turns to a massive fermion 
when ($ a ) 7^ 0. Integrating them out, will induce an 
effective non-Abelian Chern-Simons (CS) action^, i.e., 



S cs = ± J d z xe^ x [A^d v A a x + |e ai,c A»^^]. Combin- 
ing with Sr , the theory will turn into non-Abelian Chern- 
Simons-Higgs theory. The winding number of the case is 
found to be k/no and (fc + q)/no (k is an integer, uq is 
the center of the gauge group, q is a real number) for 
topological and non-topological vortices, respectively 2 ^. 
If ($ a ) = 0, the \& fermions are massless. R has a spin 
gap, so the system are in the semi-metal with a spin gap 
phase. 

The construction of non-Abelian vortices depends on 
the spin-fcrmion formula. In our analysis, we do not 
consider the physics in the valley and sublattice spaces. 
Thus, this decomposition can be applied to other types of 
lattices. Our discussion also holds for the 3 dimensional 
physics with X3 as the direction of the vortex string. 

In conclusion, we show that there is an local gauge 
redundancy in the emergent gauge theory of the Hub- 
bard model, and this gauge theory is U(2). Based on 
the effective Lagrangian in Eq.©, we construct the non- 
Abelian vortex solutions , and show that there are 2k 
orientational zero modes for fc-winding vortex. For the 
elementary vortices, there are 2k ZESs for the fermions. 
Far away from the vortex core, the system is in the in- 
sulator phase with Neel order. In the vortex core, the 
system are in the semi-metal phase with a spin gap. 
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